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ABSTRACT 

Three kinds of astrogeodetic geoid maps for Japan are presented: one referred 
to the global (18, 18) geoid of the 1973 Smithsonian Standard Earth (III) (SE HI), one 
referred to the best-fitting ellipsoid of SE ni, and one referred to the reference 
ellipsoid of the Tokyo datum. Interpolations of the deflection of the vertical are 
carried out by a least-squares estimation method. The geo id-height differences 
obtained are compared with solutions of satellite-derived station positions. Good 
agreement is found in a comparison with dopple r-tracking stations. 
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ASTROGEODETIC GEOID OF JAPAN* 

Yasuhiro Ganeko^ 

1. INTRODUCTION 

There are several ways of obtaining the geoid undulation. One is through a 
satellite-derived geoid, which is limited to long-wavelength components owing to the 
inherent inaccuracies of satellite tracking. Another is the gravimetric geoid, com- 
puted by applying Stoke's formula to surface-gravity data. This method, however, 
requires gravity values covering the whole surface of the earth, and such dense 
coverage is currently difficult to obtain. 

When a detailed and relative geoid is to be determined in a restricted area and 
when there are sufficient data in that area, the astrogeodetic geoid is practical to 
obtain, in that it does not need worldwide data. Even though observations of the 
deflection of the vertical are time-consuming, the astrogeodetic geoid is valuable 
because interpolation is possible if suitable amounts of vertical-deflection observa- 
tions exist. 


This work was supported in part by Grant NGR 09-015-002 from the National 
Aeronautics and Space Administration. 

^Visiting Scientist, Smithsonian Astrophysical Observatory, October 1974 through 
October 1975. Currently in the Hydrographic Department, Maritime Safety Agency, 
Tokyo. 
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2. AS TROGE ODE TIC GEOID 

The difference in geoid height between astrogeodetic stations is computed from 
Helmert's formula 

i+1 

N i+ ^-N. = -J* (| cos A + rj sin A) ds , (1) 

i 

where N. is the geoid height of the ith station; | and rj are the vertical-deflection 
components in the meridian and prime vertical, respectively; and A is the azimuth 
of the direction of the ith statical to the (i + l)th station. 

The current method of interpolating the deflection of the vertical needs dense 
gravity measurements around the astrogeodetic stations. Another possibility of 
interpolation, when the vertical-deflection stations are sufficiently dense, uses least- 
squares estimation (Heiskanen and Moritz, 1967). Least-squares methods also 
provide an estimate of accuracy. The number of vertical-deflection observations 
that have been made in Japan since 1886 is more than 450. This distribution is suffi- 
ciently dense for most areas in Japan. 

To find the accuracy of an estimate of the deflection, let the vector V be 
• 


Assuming that the vertical deflection at point p, V , can be estimated by the linear 
combination of observed values, V., we can write 





Ev v i - 


(2) 
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where a pi are scalar coefficients. This means we can neglect the correlation between 
| and ?]. The difference between the correct value at p and the estimated one is the 
error of estimation: 


£ = V - ‘V = V -V Q n ? 

P P P P LmJ Pi 1 
i 

By squaring, we get 




( v p-E“pi v i) ( v p 


-E“ P i \ 


= ( 5 P + ’’p) ' 2 E Q pl ( Vi + Vi* + E Q pi°pJ <Mj + Vi* • < 3 > 

i ij 

From equation (3), the root -mean-square error of the estimation, m^, is given by 
m n = + cPP - 2 V' a . ( C ?l + C pi ) + V* a. a . (C ij + C ij ) , (4) 

P 44 VV P 1 V 44 W ^ pi pj V ££ T\r\f 

i ij 

AB AB 

where C „ . and C are covariance functions of the vertical-deflection components 

44 vv 

| and tj, respectively. These covariance functions depend on distance and direction. 
Using the expression 


G, =C* B + C AB , 

ab a m ’ 


we can write 


m — G - 2 \ a . G . + N a , a . G.. . 

p pp pi pi pi p] ij 


(5) 


ij 


4 



To find the coefficients a pi that minimize the estimation error, the necessary condi 
tions are 

0m 2 

^ ~ 0 , i = 1, 2, . . . , K , 

^Pi 


or 


K 

y! a pj G ij = G pi * 1 = 1, 2, . . . , K , (6) 

j=l 


where K is the number of observations used for the estimation. Expression (6) is a 
system of linear equations with K unknowns, a p j. By use of matrix notation, equa- 
tion (6) can be written as 


^ra. “ G , 
P P 

a = G ” 1 G 
P P 


0) 


where 





r_ — 


““ ^ 


V 


g pi 


G 11 G 12 * * * °1K 


V 


G P2 


G 21 G 22 

Q = 


, G = 

• 

, G = 

• 

* ■ 

P 

• 

* P 

• 


• ♦ • 


> K _ 


• 

>_ 


• • 
g ki ■ ' ' C KK 






w 


The standard error of least-squares estimation is given by 



-( g p g " 1 g p) • 


( 8 ) 
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3. COVARIANCE FUNCTION G AB 

The observed deflection of the vertical includes both the systematic datum shift 
and the deflection included in the global geoid. To obtain the detailed local geoid by 
means of a covariance function, these two components should be eliminated. 

For tiie datum shift, we introduce the translation values Ax, Ay, and Az of the 
geodetic datum, determined by satellite tracking. According to the 1973 Smithsonian 
Standard Earth flU) (SE HI) (Gaposchkin, 1973, 1974), these values for the Tokyo 
datum arc 


Ax = -136 m , 
4 ? =+521m , 
/lz = +681 m . 


(9) 


By adding these values to the station coordinates on the Tokyo datum and neglecting 
the geoid height and passible datum rotation, the new three-dimensional geocentric 
coordinates become 


x' = (N + h) cos 4> cos X + Ax 
y' = (N + h) cos 4> sin X + Ay 

[(a) ** + k] sin 4> + Az 


z' = 


i 

t 


where 4> and X are on the Tokyo datum and h is the mean-sea-level height; N is given 
by 

n a 2 

V cos V b 2 sinV'' 2 ’ 

b-a(l-f) , 
a = 6377397.2 m , 
f= 1/299.1528 , 
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in which a and f are the parameters of the reference ellipsoid of the Tokyo datum. 
The coordinates x y', and z' can be converted to the geocentric <j/ f \' t h' system 
throug v he inverse transformation of 


x' = (N' + h') cos <J>' cos X' f 


y' « (N' + h') cos <J>' sin X' , 


2'-[(|) 2 N' + h'] e in4,' , 

b' = a'(l - f') , 

t a' 2 

(a' 2 cosV + b ' 2 sin 2 <f >') 1/2 


» 


for which the ellipsoid parameters 

a' = 6378140 m , 
f'= 1/298.256 , 


are those of the best-fitting ellipsoid of the SE HI solution. We can calculate the com- 
ponents of the deflection caused by datum shifts from the differences between the datum 
coordinates (<f>, X) and the geocentric ones (4>', X'): 


(*)_ 


w 

cos <j> (X - X') 


( 10 ) 


The components included in the global geoid are given by using the satellite- 
derived geopotential of the earth. These are 


V e 



Rg Q cos 



( 11 ) 
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where T is the satellite-derived anomalous geopotential, R is the mean radius of 
s 

the earth, and g Q is the mean value of gravity over the earth. From the SE HI geo- 
pctential coefficients C^, up through degree and order 18, plus the following 
earth's constants, 

a * 6378140 m , 
f= 1/298.256 , 

20 -2 

GM = 3. 986013 X 10 cm° sec , 
w= 7. 292115085 X 10" 5 sec" 1 , 


the anomalous geopotential, from a spherical approximation, is given by 


18 n 

T s =(t)EI < 6C nm cos + sln mX > C < sln ♦> - 

n=2 m=0 


6C = C - C , 
nm nm nm ’ 


6S = S - S = S 
nm nm nm nm 


Following the expressions of Caputo (1967), the normal geopotential coefficients 

C° , S° are 
nm* nm 


S = 0 , 
nm ’ 


nm 7 


when m * 0 or n is odd , 


C° = il±L i + 

2i,0 2i + 1 1 3(2 i 


8f K 2 1 Jt l 
+ 3)GiT *** 2_f * » 


i = i, 2, ... , 


c 3. 3 2 .. ^ c 2, 

K = £ b ( 1+ £ ) 

2 21(3 + e 2 ) tan -1 €- 3E] 


b = a(l -f) , 
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express the local variation of the geoid. These residuals were calculated for Jap-n, 

by using the S£ M geopotential for equation (11), at every station where vertical- 

deflection observations were carried out from 1947 to 1973; there were 284 such 
★ 

sites. The ir. jut; and the mean squares of the residuals are 

U R > = - o:'23 , (4> = (7:'58) 2 , 

<ri R ) = + li'64 , <4>=(8!'16) 2 . 


The next step is to evaluate the covariance function G A0 by using the vertical- 
deflection residuals. When the anomaly correlation distance is small compared to the 
radius of the earth and when the anomaly is statistically homogeneous, the following 
relation exists between the gravity-anomaly correlation and the vertical-deflection 
correlation (Shaw, Paul, and Henricson, 1969; Jordan, 1972): 


C 




+ c 

vv 




(13) 


Bulletin of the Geographical Survey Institute, Japan, vol. n, parts 2 and 3, 1951; 
vol. m, parts 2 to 4, 1953; vol. IV, parts 3 and 4, 1955; vol. V, part 4, 1958; vol. XVIII, 
part 1, 1972; and unpublished. new data. 
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where C is the covariance function of the local gravity anomaly. If the gravity 
66 


anomaly is • sotropic, then, from equation (13), G^ B is a function of distance only. 
By averaging + C in the azimuth, the covariance function G for Japan is wei: 
represented by an analytical expression of the form (see Figure 1) 


G(r) = a 2 * * * * e _r ^ , 


(14) 
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where I = 55 km and <j = (lli'l) . From equations (13) and (14), the covariance 

function of the local gravity anomaly in Japan is obtained: 

^ _ ~2 ^-r/D 

(15) 


gg' 7 g 
Or 2 = (53 mgal) 2 
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where a* is the mean-square value of the local gravity anomaly. The above analytical 
expression for C (r) gives the following covariance functions (Shaw etal. , 1969): 


C U 


c se ( r > 

— “I cos 2A f c (r) 


C = C7j 

vv d 


C„(r) 

+ cos 2A f (r) 

a 

L u g 


(16) 


c 4„ = c vr ' sin ^ f <= (r) ’ 


where A is the azimuth measured clockwise from North, 



and r is measured by the unit D. The values obtained from vertical-deflection data 
and the analytical model are shown in Figures 2a, b, and c. 

The absolute value of the computed normalized covariance of £ and rj barely 
exceeds 9. 1. Therefore, it is considered suitable to neglect the correlation between 
| and rj. 


12 



/Oi-CT, 



Figure 2c. Correlation between £ and tj. 
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4. GEOID HEIGHT 


Using Helmert's formula [equation (1)] and the least-pquares estimation method 
[equations (2) and (7)] to interpolate the deflection of the vertical, we can calculate 
the relative local geoid referred to the SE HI global geoid. The integration in equa- 
tion (1) is replaced by the summation 

^k+i + 

Xi + i-xr-22 ^-—~ s • < 17 > 

k 

r»j 

where S is chosen to be 10 km (see Figure 3), is obtained from 


cos A + sin A , 
and the estimated deflection of the vertical is obtained from equation (2). 



NORTH 



Figure 3. Geoid-height integration steps. 

The observed stations covered a square 2° X 2 ° area around the estimated point. 
The number of stations in each area varied between 20 and 40. Adding the old data 
(Torao, 1949) brings to 450 the total number of vertical-deflection stations used to 
calculate the local geoid for Japan, excluding the old data from the Hokkaido area. 
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In the Hokkaido area, data before 1947 were not used, because of inconsistencies 
between the old and the new. Relative weights were not taken into consideration 
because of the uncertain accuracy of the old data, even though the old data might 
include large errors. The purpose of this paper, however, is to examine the effi- 
ciency of the least-squares estimation method in interpolating vertical deflections. 
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5. COMPARISON OF INTERPOLATION METHODS 


The current method of interpolating vertical deflections is a gravimetric one, in 
which dense gravity observations around astrogeodetic stations are required. A test 
of gravimetric interpolation was carried out by Ono (1974) in a restricted area in Japan, 
and geoid-height (quasi-geoid-height) differences were calculated between astro- 
geodetic stations. Interpolations by means of least-squares estimation were also 
made along the same route, shown in Figure 4, and the results are presented in 
Figure 5. The geoid-height differences were transformed to those referred to the 
best-fitting ellipsoid of SE HI. The maximum difference in geoid height due to differ- 
ent interpolation methods is 44 cm, and at the farthest point, numbered 263, the 
difference is 16 cm. These results would indicate that the least-squares estimation 
method is one of the most effective ways to interpolate vertical deflections. With this 
method, we can easily calculate the relative geoid height, without additional data, 
over that area of Japan that is densely covered with astrogeodetic stations. 

140° 0 / 

+ + 4- 38°3o' 



Figure 4. Portions of the route of the two interpolation methods, astrogravimetric 

and least-squares estimation. The astrogeodetic stations are indicated by 
circles and are numbered. 


17 



QUASI-GEOID HEIGHT (m) 


5 


0-— O ASTROGRAVI METRIC 


Ar— -A LEAST-SQUARES ESTIMATION 



262 126 248 60 266 267 268 263 

36°00'N ASTROGEODETIC STAVON NUMBER 38°37'n 

140°02 / E 140° 20' E 


Figure 5. Comparison of the two interpolation methods. 


18 



6. GEOID MAPS FOR JAPAN 


We examined the closing accuracy over Japan by dividing the country into test 
loops, labeled A, B, etc. and shown in Figure 6. The absolute values of the closing 
errors are also given in Figure 6. The closing errors are ~1 m at the most and 
usually less than 1 m, and the average is 0. 7 1 m. 



Figure 6. Loops used in the testing of closing accuracy. 


To obtain the geoid undulation over Japan, series of loops of geoid profiles were 
calculated and adjusted so that the loop closures became zero. The adjustments were 
carried out by hand through trial and error until reasonable corrections could be found. 
The density of the data and an estimate of the accuracy were taken into consideration 
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in the adjustment procedure. The length of each loop was between 200 and 500 km. 
Some comparisons between adjusted and directly integrated geoid-height differences 
[by use of equation (17)] are given in Table 1, and the locations of comparison points 
and the routes of direct integration are shown in Figure 6 by broken lines. The 
shaded portions in Figure 6 indicate where large corrections had to be added, in 
mountainous and s parse-data regions. More vertical-deflection observations are 
needed in these areas to achieve the same accuracy obtained from geoid-height 
calculations. 


Table 1. Comparison of adjusted and nonadjusted geoid-height differences. 


Beginning point 
of integration 

End point of 
integration 

Direct 

integration 

(m) 

After 

net adjustment 
(m) 

Direct minus 
adjusted 
(m) 

P 

Q 




4.= 33°30'N 

35°00’ 

-1.98 

-0.96 

-1.02 

\= 131°30'E 

138°00' 




Q 

R 




4> = 35°00' 

38°00' 




X = 138°00' 

141° 00' 

+3.05 

+2. 87 

+0. 18 


The local geoid for Japan, calculated from the reduced vertical deflections 
and rfo and referred to the SE m global geoid reproduced in Figure 7, is shown 
in Figure 8. The geoid height at the datum point is taken to be zero. The contour 
shapes are quite similar to those of the gravimetric geoid by Hagiwara (1967) (see 
Figure 9), although there are large discrepancies in the contour values themselves. 
These discrepancies, however, seem reasonable, given the fact that the gravimetric 
geoid was calculated from gravity data in a restricted area. Hagiwara's gravimetric 
geoid may be considered representative of local variations of the geoid, perhaps 
including large truncation errors (Hagiwara, 1970). 
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120° E 130 140 150 160 

Figure 7. SE HI geoid through degree and order 18, f - 1/298.256 (in meters). 

Adding the global geoid (Figure 7) to the local geoid (Figure 8), we have the 
geoid referred to the SE m best-fitting ellipsoid; this is reproduced in Figure 10. 

The geoid height at the datum point is again taken to be zero. The geoid referred to 
the reference ellipsoid of the Tokyo datum can also be obtained through transforma- 
tion by using translation values from equation (9) and ellipsoid parameters. This 
type of geoid is shown in Figure 11. Fisher (1960) published the same kind of astro- 
geodetic geoid for all of Japan except for Hokkaido and the southern part of Kyushu. 
The geoid-height differences between hers and Figure 11 are as much as 3 m in some 
areas. 
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7. COMPARISON OF STATION COORDINATES 


The geoid-height differences referred to the SE HI best-fitting ellipsoid, 
a = 6378140 m, f = 1/298.256, were compared with solutions derived from SE m and 
from doppler observations (Anderle, 1974). The results are presented in Table 2. 
Locations of the satellite-tracking stations are shown in Figure 12. 

Table 2. Comparison with satellite-derived geoid heights. 


Geoid-height differences 


Stations 

(geoid height, m) 

Distance 

(km) 

Astro- 

geodetic 

(m) 

SE m 
(m) 

* 

Doppler 

(m) 

Height 

discrepancy 

(m) 

Tokyo - Dodaira 
(35.2) (38.3) 

48 

-3.1 

-3. 1 

(constrained) 

- 

0.0 

Tokyo - Kanoya 
(35.2) (22.8) 

933 

-6.0 

-12.4 

(partly 

constrained) 

- 

6.4 

Sasebo - Misawa 
(30.3) (31.7) 

1334 

-2.3 


-1.4 

0.9 


* 

Anderle (1974; private communication, 1975). 


The geoid height of Kanoya, which is 22. 8 m in the SE HI solution, is somewhat 
too low, considering the local geoid of Japan. The large discrepancy between the 
SE HI and the astrogeodetic-geoid solutions seems to arise partly from the adopted 
geoid height (-19 m) of Kanoya on the Tokyo datum in SE El, in contrast to the -12 m 
determined in the astrogeodetic geoid. 
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Figure 12. Locations of satellite-tracking stations. The satellite-derived geoid 
heights (in meters) are shown below the station names. 




8. CONCLUSIONS 


The efficiency of the least-squares estimation method in interpolating deflections 
of the vertical has been shown. This method requires only vertical-deflection data. 
When these data are dense enough - roughly speaking, more than 30 in a 2° X 2° 
area - this method is one of the most effective ways of calculating geoid undulation. 

The correlation distance of the local anomaly field in Japan is 55 km, and the 
root-mean-square value of the deflection of the vertical of the local anomaly field 
is lll'l. 

Good agreement is found in the comparison of doppler results with the geoid 
obtained here. 
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